Nanomechanics of an individual vortex in an anisotropic type-II superconductor 
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As shown in recent experiments [Auslaender et al, Nature Physics 5, 35 (2009)] magnetic force 
microscopy permits one not only to image but also to manipulate an individual vortex in type-II 
superconductors, and this manipulation provides a new powerful tool to study vortex dynamics and 
pinning. We derive equations that describe the deformation of an individual vortex in an anisotropic 
biaxial type-II superconductor under the action of the microscope's magnetic tip. These equations 
take into account the driving force generated by the tip, the elastic force caused by the vortex 
deformation, and the pinning force exerted by point defects. Using these equations, we reproduce 
the main features of the experimental data obtained by Auslaender et al. 

PACS numbers: 74.25.Qt, 74.25.Sv 



I. INTRODUCTION 

In a recent paper— magnetic force microscopy (MFM) 
was employed to image and manipulate individual vor- 
tices in a single crystal YP^CuaOg.ggi, directly mea- 
suring the interaction of a moving vortex with the local 
disorder potential. Several unexpected results were ob- 
tained in that paper. In particular, the authors of Ref. Is 
found a dramatic enhancement of the response of a vor- 
tex to pulling when they wiggled it transversely. In addi- 
tion, they discovered enhanced vortex pinning anisotropy 
in this crystal. These results demonstrate the power of 
MFM to probe microscopic defects that cause pinning 
and show that the described manipulations of an indi- 
vidual vortex provide a new powerful tool for studying 
the vortex dynamics and vortex pinning in type-II super- 
conductors. 

In this paper we derive equations that govern the vor- 
tex dynamics under such MFM manipulations, and by 
solving these equations numerically, we provide some in- 
sight into the results of Ref. QJ 



II. EQUATIONS FOR A MOVING VORTEX 

Consider a platelet-shaped biaxial anisotropic super- 
conductor, with its crystalline c-axis being perpendicu- 
lar to the plane of the platelet (and the a and b axes 
in this plane). Let there be a vortex directed along the 
c-axis in the sample. We denote this axis as the z axis, 
and choose the x and y axes along the a and b axes of 
the crystal. MFM employs a sharp magnetic tip placed 
near the surface of the platelet. The tip magnetization 
exerts an attractive force F on the vortex end. This force 
can shift the top of the vortex when the tip moves. On 
the other hand, it is possible to measure dF z /dz at the 
tip, and this permits one to visualize the position of the 
top end of the vortex^ Let X, Y be the position of the 
tip in the x-y plane, while its height above the surface 



of the platelet be Z. We shall describe the shape of the 
vortex by the functions x{z) and y(z) with z < 0, the po- 
sition of the vortex end at the surface is thus xq = x(0), 
y = y(0). Below we shall use the following dependence 
of the force F on height Z and on the two-dimensional 
vector R = (A" — x , Y — y ): 
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F = q 



R+(Z + fe )z 

(r 2 + {z + M 2 ) 3 / 2 ' 



(1) 



where the constant ho f» A (A is of the order of the Lon- 
don penetration depth), q = m$o/27r, <!>o is the flux 
quantum, m is the magnetic monopole strength of the 
tip (or the magnetic moment per unit length of a long 
narrow cylinder used as tip), and z is the unit vector 
along the z axis. This dependence is obtained if one con- 
siders the tip and the end of a straight vortex as magnetic 
monopoles of strengths rh and 2<&o / /^o-— The lateral com- 
ponent of F, Fi at , gives the driving force acting on the 
vortex. The maximum of Fi at with respect to variations 
of R is reached at R = [Z + ho)/y/2 and is equal tc£ 
F m w 0.38bq/ (Z + ho) 2 . In our following numerical cal- 
culations we shall use formula (|TJ) even when the vortex 

is curved, and to describe the lateral component fexdz 
of the external driving force applied to a vortex segment 
which has the projection dz on the z-axis, we shall em- 
ploy the model expression 
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(R 2 + (Z + ho) 2 ) 3 / 2 A 



(2) 



This expression can be justified if the change of the total 
lateral force F/ at on the scale A in the x-y plane is rela- 
tively small (i.e., if R ^> A). However, in the case when 
the vortex shift (x§ + yo) 1 ^ 2 caused by the tip is essen- 
tially larger than A, this shift is practically independent 
of the specific form of the z-dependence of f]L, see below. 
So, to clarify the physics without additional mathemat- 
ical complications, below we shall always use the model 
dependences (J]) and ([2]). 
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As it was mentioned above, measurement of dF z /dz is 
employed to visualize the vortex. Equation fl} yields the 
following expression for this derivative: 



dF z 



dz 



\2(Z + h ) 2 



R 2 



(R 2 + {Z + h ) 2 ) 5 / 2 ' 



(3) 



This derivative is maximum \dF z /<9z| max = 1qj(Z + /irj) 3 
when the tip is just above the vortex, i.e., when X = xo, 
and Y — yQ. On the other hand, the maximum lat- 
eral force occurs when R = (Z + ho)/V%, and hence 
\8F z /dz\ = q(2/3) 3 / 2 /(Z + h ) 3 » 0.27\dF z /dz\ max at 
this R. In other words, the maximum of the lateral force 
and the maximum of \dF z /dz\ occur at different positions 
of the tip and of the vortex end. 

We shall consider the vortex as an elastic string. In 
the case of a biaxial superconductor the line tension of 
the vortex, £i{6, if, tp), and the pinning force acting on its 
unit length, f p (9,if,ip), were calculated in Ref . |6|. The 
angles 9 and f define the direction of the vortex, i.e., we 
shall describe this direction by the unit vector 



(sin 9 cos if, sin 9 sin f, cos ( 



(x',y',l) 



v 



(4) 



while the angle tp defines the direction of the pinning force 
or of the vortex distortion in the plane perpendicular to 
the vortex, Fig. 1. Here the prime means d/dz. In the 
subsequent analysis the line tension will be required only 
for the case 9 = since the linear elasticity theory is 
valid up to sufficiently large angles 9 if the parameter e 
is small X Then, we havei^ 



s t (if, tp) = e e 2 r](if + ip)=si(if + tp), 



(5) 



where e = A a b/A c is the parameter of the anisotropy; 
£o = ($o/A ah ) 2 ln(A a6 /£ af) )/(47i7io); Kb = VKM; A c , A q 
and Af, are the London penetration depths, C — A a /A& is 
the parameter of the anisotropy in the a-b plane, and 



r)(ip) = (cos 2 If + C 



_1 sin 2 if. 



(G) 



Since at 9 = the plane perpendicular to the vortex 
coincides with the x-y plane, the combination ip + tp in 
Eq. is the angle defining the direction of the vortex 
distortion in this plane relative to the x axis4 As to the 
pinning force, it is described by the expression 6 



f P {e,<p,i>) = f. 



£ab COS ( 



p ao,f^y 



(7) 



where /° is the pinning force for the vortex along the c 
axis in the uniaxial superconductor with the same A a (, 
and £ a (, = y/£ a £b- Here £ Q and £b are the coherence 
lengths, and 

£, 2 (9,f,tp) = £ 2 b (,(sinf cos ip cos 9 + cos p siml>) 2 (8) 
+ i(cos<y9Cos V>cos# — sin 93 sin -0) 2 
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FIG. 1: Definition of the angles 6, ip and tp. The angles 9 and 
ip specify the direction of the vortex shown as bold solid line. 
The angle tp in the plane perpendicular to the vortex defines 
the direction of the pinning force; tp is measured from the line 
that is the intersection of this plane with the plane containing 
the vortex and the z-axis. 



In YBa 2 Cu 3 6 99 one has£ e m 1/7 (i.e., e 2 < 1) and±>2 
Cwl.3. 

Consider a vortex segment limited by the planes z and 
z + dz and specified by the angles 9 and tp. Let us analyze 
the balance of the driving, the pinning, and the elastic 
forces applied to this segment. All these forces are per- 
pendicular to it. However, to find the two functions x{z) 
and y(z) describing the vortex, it is convenient to carry 
out the analysis in the x-y plane, projecting all the forces 
onto this plane. The projection of the elastic force acting 
on this segment, fj^dz, can be described by the simple ex- 
pression fgfdz = {ei x x" , £i y y")dz even at sufficiently large 
9 since the linear elasticity theory is valid up to the angles 
satisfying e 2 tan 2 9 <C lX Here ei x = e/(0) = £n£ 2 C and 
£/y = £/(f/2) = £n£ 2 /C are the appropriate line tensions 
at 9 = 0, see Eqs. (JS|) and ©. Adding this projection of 
the elastic force to the external force defined by Eq. J2]), 
one obtains the projection f^dz of the resultant force fdz 
on the x-y plane. Then, the first of two equations for 
x(z) and y(z) is 



f(x,y,X,Y) = f\\ 



(9) 



where /i is the absolute value of the projection of the 
so-called critical forced on the x-y plane. This critical 
force is the force at which the vortex starts to move. It is 
determined by the pinning force, but in the anisotropic 
superconductor it can differ from the pinning force. 6 Note 
that we write one equation which connects the abso- 
lute values of f " and fj rather than two equations for the 
x and y components of these forces. This is due to the 
fact that the direction of the pinning force (and hence of 
the critical force) is not known in advance and is dictated 
by the direction of f " . 

The critical force /J is determined by the following 
formulas: Let the direction of the force f be specified 
by the angle ip in the plane perpendicular to the vortex. 
This angle can be expressed in terms of the component 
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fx and fy of the force f " as follows: 



tan?/; 



cosfl(/J - fl t&nip) 
/J+/Jtan<p 



(10) 



The pinning force f p in the direction tp is given by 
Eqs. (J7J) and ©, while the critical force f c in this di- 
rection ip is determined byS 



tan(?/> — ipi) 



.WOT 



(ii) 

(12) 



where the prime means d/dtpi, and the angle ip\ in the 
plane perpendicular to the vortex defines the direction 
along which the vortex starts to move when the force 
acting along ?/> exceeds f c . The fact that ipi in general 
differs from tf> is due to the anisotropy of the pinning. On 
determining ipi from Eq. (|lip . one then finds f c (jp) from 
formula The explicit form of Eqs. (fTTJ) and (JT^Jl for 
the case of the pinning force described by formulas ([7]) 
and (JSJ) is presented in Appendix O Finally, when the 

critical force f c (ip) is found, its projection /J on the z-y 
plane is determined by the formula 



(cos 2 cos 2 ?/> + sin 2 ip) 1 / 2 



cos 6* 



(13) 



that follows from geometrical considerations. 

Equation (J5|) is a differential equation since it contains 
the derivatives x"(z) and y"(z) originating from the elas- 
tic force. As in Ref. [l], we shall consider only sufficiently 
thick superconducting crystals in which the vortex as a 
whole does not shift, and only its upper part (0 > z > zo) 
adjoining the x-y surface moves, while the lower part 
(z < zo) is pinned. The boundary point zq of this upper 
part is determined by 



x ( z o) = y(zo) = 0. 

Then, the boundary conditions to Eq. ([9]) are 



x'(z ) 

AO) 



y'(zo)=o, 

y'(Q) = 0. 



(14) 



(15) 
(16) 



If these conditions were not fulfilled, the derivatives x' 
and y' would be discontinuous at the points z = zq 
and 2 = 0, and the elastic force (si x x", £i y y") would 
be singular there J£ In the most interesting case when 
x o + Vo ^ ^ 2 ( an( l hence \z \ ^> A), one can put A — > 0. 
In this limiting case the driving force F is applied to the 
vortex only at its surface point (xo,yo)- Then, in equa- 
tion (JHJ) the force f| x can be omitted, f " coincides with 
f|, and the driving force F only modifies the boundary 
condition (TIT))) . Now the integration of the forces over the 
thickness of the surface layer gives 



Ao) = —, y'(o) = ^- 

Six Sly 



(17) 



This result shows that at small A the vortex dynamics is 
practically independent of the distribution of the driving 
force F over the surface layer of thickness A. 

Equation alone is not sufficient to find the two func- 
tions x(z) and y{z). We now derive a second equation for 
these functions. When the position of the tip changes, 
the vortex begins to move in the direction ipi in the plane 
perpendicular to the vortex. This movement of the vor- 
tex in the perpendicular plane corresponds to its shift at 
an angle ipi (measured from the x axis) in the x-y plane. 
A geometrical consideration shows that this angle ipi is 
determined by 



tan?/>i 



tan ip + cos 9 tan ip\ 
1 — cos 6 tan tp tan ipi 



(18) 



Thus, changes of the functions x{z) and y(z) in time are 
connected by the relation 



dy ~ dx 



(19) 



This is a second equation for the functions x(z) and y(z). 
Since the time t can be expressed in terms of the known 
functions X(t), Y(t) that describe the shift of the tip, 
Eq. (Til?)) and its solution (i.e., the shape of the vortex 
at some moment to) depend on the trajectory Y(X) of 
the tip in the x-y plane at previous times (t < to) rather 
than on a specific form of the temporal dependences X(t) 
and Y(t). This situation is reminiscent of the case that 
occurs in the theory of the critical states of type-II su- 
perconductors when the external magnetic field H a ap- 
plied to a superconducting sample changes in a complex 
way . 11 ! 12 In this case the critical states are different for 
different histories H a (t) with the same final value of H a . 

Equations (TTjl- lfT^)) describe the vortex dynamics in 
thick superconducting crystals when the tip moves in its 
x-y plane. We solve these equations in the next section. 



III. RESULTS 

The equations of the previous section show that if the 
driving-force density f\ x at the surface of the supercon- 
ductor, z = 0, is lower than a certain threshold /1(a) 

where a is the angle of f]L relative to the x axis, the vor- 
tex remains pinned, i.e., x{z) = y(z) = 0. In particular, if 
the driving force acts along the x or y direction, we obtain 
the following thresholds: /^VC an d fp/VC: respectively, 
which coincide with the appropriate pinning forces. Here 
we have used the formulas of Appendix A and the fact 
that S = l/( 2 > 1/2 at the experimental valued of 
( = 1.3. Equivalently, these threshold conditions can be 
rewritten in terms of the total forces, F x < F px = fp\\/C, 
and F y < F py = /pA/^/C- If the driving force exceeds the 
threshold values only a little, i.e., if F x — F px <C F x , or 
F y — F py C -R,, we find from the equations that xo or yo 
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begins to deviate gradually from zero, 



x 



2\{F X - F. 



yo 



2X(F y 



F ) 3 



SlyF* 



(20) 



With further increase of the driving force, at F x ^> F px 
or F y 3> F py but at the same time under the conditions 
F x < six = Ce 2 £o or F y < ei y = e 2 e /(, we arrive at 



x 



Fx {Fx ^-Fp X ) 

2C 3/2 / P c 



e 2 e 



yo 



S r y\ r y ^ 1 p: 
2/ p c £ 2 e 



(21) 



The additional conditions F x <C £i x , F y <C £i y mean that 
the characteristic tilt angle 9 of the vortex is small [see 
Eqs. (fl~7|) in which x'(0), y'(0) are just equal to tan(9]. 
This smallness of 9 was assumed in Ref. ll| in analyzing 
the vortex dynamics, and formulas (|21|) coincide with 
those obtained in that paper. However, F px /£i x , Fpy/^iy 
are not necessarily small in an experiment. In this case 
formulas ([2Tj) , strictly speaking, have no region of appli- 
cability. Moreover, boundary conditions (fl~7|) show that 
the characteristic tilt angle 9 is not small at typical exper- 
imental values of F x . y ~ 5 — 20 pN even when ei xy ~ 10 
pN. So we do not assume in this paper that tan0 <C 1. 
The equations of the previous section have been derived 
only under a weaker condition e 2 tan 2 9 -C 1. But when 
6 *~ 1, the critical force / c differs from f p even for symme- 
try directions.— For example, when the tip moves along 
x and thus the vortex also bends along this direction, 
formula (|A5[) of Appendix A gives 



(22) 



/J(0)=/pVC, tan 2 0<^-l; 



f!(0) = 



£3/2 



COS 9 y/C 2 - cos 2 9, tan 2 6>>— -1, 



c 2 



while for the tip moving along the y axis, one has 
/J(0) = ^|, tan 2 0<2C 2 -l; 



(23) 



/J (61) = 2C 1/2 /; cos Oy/l - C 2 cos 2 0, tan 2 6> > 2( 2 - 1. 

In other words, even at moderate 9 the critical force be- 
gins to depend on this angle, and the formulas for xq and 
y Q become more complicated than Eqs. (I21[) in which /J 
was assumed to be constant and to coincide with the 
appropriate pinning force, /J (9) = f£\/C at ^ = and 

fc(0) = f^/VC at tp = tt/2. Such a dependence of /J(0) 
in general causes that the ratio yo/%o at large driving 
forces differs from the value £ 3 ~ 2.2 that follows from 
formulas (f2"Tj) . This may lead to an imitation of the 
enhanced pinning anisotropy observed by Auslaender et 
al^- see below. 

During its motion the vortex lags behind the mov- 
ing tip until the maximum lateral force is reached at 
r m = max(iQ +y 2 ) x / 2 . At small driving force the vortex 
will remain at this r m , whereas at large driving forces 
the vortex will, in fact, partially recede after the tip has 




FIG. 2: Dependence of the maximum shift of the vortex end, 
r m = max(io + Vo) 1 ^ 2 , on the driving force F m x or F my when 
the tip moves along the x axis (ip = 0) or the y axis (tp = 
7r/2). The forces are measured in units of the line tension 
Ei X y = £ 2 £o, the lengths in units of A, f = 1.3, and P = 
fpX/ei X y — 0.5. The dashed lines show the appropriate X 
and Y positions of the tip, R m = {X 2 + Y 2 ) 1/2 , at which 
the derivative dF z /dz reaches its maximum on the returning 
path. As an example, the V-dependence of this derivative 
at (j) — 7r/2 and F my /ei xy ~ 2.2 is presented in the inset. 
The circles in the inset mark the virgin curve, and the arrows 
indicate the direction of the tip motion. 



moved away. Experimentally the final location of the 
vortex is evaluated on the returning path of the tip by 
monitoring the tip location R m = (X 2 + y 2 ) 1 / 2 at which 
dF z /dz is maximum when the tip is above the vortex (or 
closest to it). In Fig. [5] we show the maximum shift of 
the vortex end, r m , in the forward direction and R m on 
the returning path of the tip vs. the driving force when 
the tip moves either along the x axis or along the y axis. 
In these cases the vortex shifts along these symmetric di- 
rections, too. Figure [2] shows that at low driving forces 
the experimentally determined R m accurately reproduces 
the maximum shift of the vortex r m whereas at higher 
forces R m slightly underestimates the actual r m . 

In the construction of Fig. [21 as well as Figs. [3] and dj 
we put C = 1.3 and measure forces in units of the line 
tension ei xy = (sixSiy) 1 ^ 2 = £ 2 £o, and lengths in units 
of A (hence the force densities f p , f c , f e i, and fj} x are 
in units of ei X y/X). Then, taking into account the model 

dependence for the driving force density j\ x , one finds 
that equations (O and (fl~9)) for x(z) and y(z), as well as 
the boundary condition (fT7|) . become independent of the 
absolute values of Eq, e, and A. They depend only on 
the dimensionless forces F Xt y/e 2 Eo and the dimensionless 
parameter P = f^\/e 2 £Q. Thus, in a certain sense Fig. [2] 
is universal. But in this scaling procedure one has to bear 
in mind that if one changes the parameter A keeping a 
fixed value of F X!V /e 2 eo, this leads to a change of the 
tip position X, Y which is not scaled with A, see Eq. ([TJ. 
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FIG. 3: Polar plot of the maximum shift r m of the vortex 
end r = (xq + yo) 1 ^ 2 versus the angle ip of the tip motion 
(solid line). The tip moves along a straight line in the X-Y 
plane with a sufficiently large amplitude and at the height Z 
that leads to F m /ei xy ~ 2.2. The dashed line shows the shift 
R m = (X 2 + Y 2 ) 1/2 of the tip when the derivative dF z /dz 
reaches its maximum. Length unit is A, P — 0.5, £ = 1.3. 

However, if one is interested only in the tip position when 
it is just above the vortex, the scaling still holds in this 
case. On the other hand, when the relative positions 
of the tip and of the vortex are essential (e.g., in the 
construction of Figs. 151 -IT2]) . we use the following set of 
input parameters: 

A=0.2/im, Ei xv ={e lx eiy) 1,2 = 9pN, 

P=^—=0,5, — =l.l/im 2 , Z+/i =0.44/im. (24) 

The data of Fig. [5] are similar to the data of Fig. 3b 
in Ref. [l|. Moreover, a semiquantitative agreement of 
these data can be obtained if one takes A of the or- 
der of several tenths of a micron and Ei xy ~ 10 pN. 
However, this value of the line tension ei xy is 10 — 20 
times larger than the theoretical estimate of this quan- 
tity, e ixy = £ 2 ($o/A a6 ) 2 ln(A Qb /£a b )/(47r/! ), at e = 1/7, 
A a 5 = 0.2 fixa, m(A a b/£ a b) = 4. Thus, apart from an en- 
hanced anisotropy of pinning discovered by Auslaender et 
oi.ji their experimental data in fact means that either the 
vortex has an enhanced line tension, or the model depen- 
dences Q and ||5J) for the driving force are oversimplified 
under the conditions of the experiment and lead to an 
essential overestimation of this force. 

In Fig. [3J that is similar to Fig. 4c of Ref. [l], we 
show the dependence of the maximum shift of the vortex 
end, r m , on the angle ip at which the tip moves along 
a straight line in the x-y plane at a certain height Z 
above the surface of the superconductor. This height 
determines the maximum driving force F m applied to 
the vortex, and in Fig. [3] this height is chosen so that 
Fm/zixy ~ 2.2. For comparison, we again show the po- 



FIG. 4: Path of the vortex end x(t), y(t) when the tip os- 
cillates along the diagonal X(t) — Y(t) (dashed line) with a 
large amplitude; same data as in Fig. [3] Both the tip and the 
vortex start at x = y = 0. Length unit is A, P = 0.5. The 
vortex path cycles a narrow hysteresis loop as indicated by 
the arrows. Due to the in-plane anisotropy £ = 1.3, this loop 
is tilted away from the tip-path (</? = 7r/4) towards the y axis. 
Also depicted is the six times enlarged and shifted path near 
the first and the second turns. The dots on the curves are at 
equidistant times. 

sitions of the tip, R m = (X 2 + Y 2 ) 1 / 2 , at which the 
derivative dF z /dz reaches its maximum. The anisotropy 
of the vortex shift, r m (tp — w/2)/r m (ip = 0) « 2.5, 
seen in the figure approximately coincides with the ra- 
tio Rm((p = n/2)/R m (<p — 0), and at £ = 1.3 this 
anisotropy is lower than the appropriate experimental 
value ~ 3.5ii This experimental value can be fitted if 
one takes C = 1-43. Thus, although this £ = 1.43 ob- 
tained with taking into account the ^-dependence of /J 
is less than £ = 1.6 derived in the simplified analysis 1 , our 
approach still cannot completely describe the enhanced 
anisotropy of pinning within the framework of collective 
pinning theory by point defects. Auslaender et o/.— sug- 
gested that this enhanced anisotropy is due to a cluster- 
ing of the point defects. 

Interestingly, when the tip moves along a straight line 
different from the x and y axes, the trajectory of the 
vortex end performs a "hysteresis loop" with its axis de- 
viating from the direction of tip motion, Fig. [4] Also 
depicted in Fig. [4] is the six times enlarged path near 
the first and the second turns, showing that the vortex 
end reaches maximum elongation, then it recedes when 
the tip moves away, and when the tip returns, the vortex 
end approaches the tip and reaches maximum elongation 
a second time. These results clearly demonstrate that the 
vortex in general moves in a direction different from the 
direction of the tip motion, and that the vortex position 
depends on the trajectory of the tip at previous times. 

In Ref. [l] the derivative (dF z /dz) was measured when 
the tip oscillates with a large amplitude along some line 
and at the same time it is slowly shifted in the perpen- 
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FIG. 5: The zigzag path xo(i), yo{t) of the vortex end (bold 
lines in the center) when the tip oscillates with large ampli- 
tude a = 1.6 /im along x and at the same time moves slowly 
up along y, with = 1/80. Tip and vortex start at 

x = y = 0. Length unit is /j,m, £ = 1.3, the other parame- 
ters are listed in Eqs. (I24|l . The aspect ratio of this path is 
max(i/o)/max(xo) ~ 3.7. The almost horizontal dotted lines 
at equidistant y = yi show the tip path when it moves from 
the left to the right (see arrows) and serve as zero lines for 
the force derivative g(x,yi) = dF z /dz plotted versus x as 
yi + 0.2 • G(x,yi) (solid lines) with G = g/max(\g\). Note 
that these curves are asymmetric due to the unidirectional 
tip motion shown here. The return path yields similar curves, 
obtained from the depicted curves by the reflection x — > —x. 



dicular direction. In this case an enhanced shift of the 
vortex along the slow scan direction was discovered, see 
Figs. 1 and 2 in Ref. [J We have investigated this situa- 
tion theoretically. In Fig. [5] the zigzag path Xo(t), yo(i) 
of the vortex end is presented when the tip oscillates 
with a large amplitude along x and at the same time 
moves slowly up along y. We also show the X profiles of 
(dF z /dz) at various fixed values of Y . Note that these 
profiles are asymmetric and are different for tip motion 
from left to right and from right to left. The data of 
Fig. [5] qualitatively reproduces the experimental data,* 
Interestingly, this figure also clearly shows how the elas- 
tic force drags the vortex back towards the origin when 
the tip goes far away from the vortex. 

In Fig. [6] we compare the vortex paths for various ratios 
of the scan rates along x and y. It is seen that the results 
are different for different rates although we do not take 
into account the effect of vortex creep here. This differ- 
ence in the vortex paths is due to the above-mentioned 
dependence of the vortex position on the trajectory of the 
tip at previous times. In this figure we also present the 
vortex-shape functions x(z) and y(z) at some moment of 
time. These functions show that during the zigzag mo- 
tion the vortex is bent and twisted into a complicated 
shape. The lower part z < zq of the vortex is rigidly 
pinned (has exactly x = y = 0) and at the surface z = 
the vortex ends perpendicularly. We find that for the tip 
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FIG. 6: The zigzag path xo(t), yo(t) of the vortex 
end as in Fig. [5] but at Y/\X\ = 1/40 (left plot, 
max(yo)/nrax(a;o) ~ 3.6) and at = 1/160 (middle 

plot, max(yo)/max(xo)=4.1). The right plot shows the vortex 
shape expressed as x(z) (solid line with circles) and y(z) (solid 
line with dots) at the moment when xo = 0.05, yo = 0.3 in 
the left plot. The dashed lines show these functions at three 
previous time steps. 
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FIG. 7: The zigzag path xo(t), yo(t) of the vortex end as in 
Fig. [S] but for C = 1 (left plot) and £ = 1.5 (right plot). The 
aspect ratio max(i/o)/max(a;o) is approximately 2.2 for f = 1 
and 5.5 for £ = 1.5. 



motion of Fig.O at z > zq the component y(z) increases 
with z monotonically, while x(z) after several scan peri- 
ods exhibits strongly damped oscillations. 

In Fig. [7] we analyze the dependence of the zigzag vor- 
tex motion on the anisotropy parameter £. It is clear from 
the figure that the shift of the vortex end in the slow 
scan direction and the aspect ratio max(j/ )/uiax(a;o) 
increase^ with increasing £. But importantly, even in 
the case of isotropic pinning in the x-y plane, i.e., at 
C = 1) this aspect ratio remains considerably larger than 
unity. From a qualitative point of view, this enhanced tilt 
of the vortex along y is caused by the fact that during the 
zigzag motion the vortex predominately moves in the x 
direction, the pinning force is also directed mainly along 



7 



o 









JT ' " ' ' ' 3> *** 








^> " ' parabola j*^' 




x(z) 












, ' ' parabola 






y(z) 






f straight .gP" 


"1.5 




















">. 


y'(z) / / \ 








. 1 




J 


■' '" sT 


"x 












/ f x'(z) I 








0.5 






: z o 





-0.6 -0.4 z -0.2 

1 1 1 1 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 

x, y 



FIG. 8: The vortex shape during oscillations of the magnetic 
tip above a superconductor that is isotropic in the a-b plane 
(£ = 1). Similar case as the left plot of Fig. [7] but to clar- 
ify the situation, we take X a b = 0.05 (im and P — 0.25 here. 
Shown are the maximum vortex displacement x(z) at the first 
excursion of the tip (i.e., at Y = 0) and the maximum dis- 
placement y(z) at the moment when xo = while t/o reaches 
its maximum value after many tip oscillations. The dash- 
dotted straight line reveals that the curve y(z) has a long 
zero-curvature segment, see also x'(z) and y'(z) shown in the 
inset (the small hump seen in the flat part of y'{z) oscillates 
in time). At small x and y both x(z) and y(z) wee parabolas 
with curvature fp/eixy (dashed lines). 



x, and hence this force opposes only the vortex tilt in the 
x direction. These considerations are supported by the 
data of Fig. [5] in which for the case of a small A we show 
x(z), the maximum displacement of the vortex when the 
tip moves only along the x axis (i.e., during the first oscil- 
lation of the tip in the left plot of Fig. [7J , and y(x) at the 
moment when j/o reaches its maximum value after many 
oscillations of the tip. Since at small A the driving force 
concentrates near the surface of the superconductor, in 
the bulk of the sample the elastic force associated with 
the curvature of the vortex has to be balanced mainly by 
the pinning force. Then, the long straight segment of the 
line y(z) shown in Fig. [5] means that the y component of 
the pinning force is practically absent in this segment. 

Some insight into the origin of the vortex-motion 
anisotropy seen in Fig. [7J can be also obtained from a 
simplified two-dimensional (2D) model. In this model 
we disregard the dynamics of the entire vortex and 
consider only the vortex end as a point (xo,yo) elasti- 
cally connected to the origin of the x-y plane, F e i — 
— (k x xo, k y yo), where k x and k y = fc x /£ 2 are some spring 
constants modelling the elasticity of the vortex. In this 
simplified approach the problem of the vortex motion be- 
comes two-dimensional, and instead of the force densities 
we deal with the elastic force F e j , the total pinning force 
F p , and the driving force Fi at — (F x ,F y ) determined by 
Eq. UJ) . The balance of these three forces and the vortex- 
end motion can be still described by the equations of 
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FIG. 9: Path of the vortex end and force balance for a sim- 
plified 2D model, see text. The tip moves as in Fig. [7J the 
left plot. Here Z + h = 0.44 fim, q = 9.9 /im 2 -pN (which 
gives F m w 20 pN); F p = F m /4; C = 1; = k = 32 
pN//^m; x and y are measured in um. The aspect ratio 
r = max(j/o)/max(yo) ~ 1.24. The force balance is shown 
for the point (xo,yo) = (0, max(yo))- 

Sec. II if one puts 8 = (p = and replaces the force den- 
sities by the total forces in the equations. Interestingly, 
in this simplified 2D approach one can qualitatively re- 
produce the main results which have been obtained above 
accounting for the real 3D shape of the vortex. In par- 
ticular, in Fig. [9] we show the zigzag path of the vortex 
end in the case £ = 1 (isotropic elasticity and pinning in 
the x-y plane). In the construction of this figure we use 
the same parameters for the tip as in Fig. [7] (i.e., we have 
F m w 20 pN). Besides this, we take F p = F m /A m 5 pN. 
This relation also corresponds to the case of Fig. [7] if one 
implies that F p for the 2D model is equal to f°\. Such 
choice of F p is dictated by a comparison of the conditions 
F m > F p and F m > f p X for a vortex to start to move 
in the simplified 2D model and in the three-dimensional 
theory. The spring constant in Fig. [5] is chosen such that 
max(xo) is the same as in the left plot of Fig. [7] The vor- 
tex trajectory presented in Fig. [9] reveals the anisotropy 
of the vortex motion in the y and x directions with the as- 
pect ratio r = max(j/o)/max(so) ss 1.24. This anisotropy 
can be understood from the following simple consider- 
ations: The maximum displacement of the vortex end 
along x is found from 

max(x ) = Fm - Fp , (25) 

where F m is the maximum value of the driving force and 
k = k x . The yo reaches its maximum when xq ~ 0, 
the vortex-end velocity v is practically parallel to x, and 
thus the pinning force is along this axis, too, Fig. [9l The 
driving force at this moment is maximum, F — F m , and 
is directed at an angle a with respect to the x axis, while 
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FIG. 10: Path of the vortex end when the tip oscillates along 
x with large amplitude and moves down from large positive 
Y > 2 to large negative Y < —2 (left plot) and then moves 
up again to large positive Y > 2 (right plot). The straight 
vortex waits at x — y — 0. When the tip approaches from 
above, the vortex end suddenly jumps to the tip and starts to 
oscillate with large amplitude, following the tip downwards. 
After some time the vortex end comes to a halt as in Figs. [5] 
|6l and [7] When the oscillating tip approaches again from 
below, the vortex end starts to oscillate with slowly increasing 
amplitude along a path that looks similar to the path on which 
the vortex end came to a halt. The vortex paths shown at 
the lower left and at the upper right are nearly identical. The 
parameters are as in the left plot of Fig. [6] 
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FIG. 11: Path of the vortex end for parameters as in the left 
plots of Figs.|6]and 1101 but now the tip oscillates along y and 
approaches the vortex end (that waits at x — y = 0) along x 
from far left (X < — 2), moving further until the vortex end 
comes to a halt. 



the elastic force acts towards the origin. Then, the force 
balance for the x and y components gives 



Fm cos a 



F m sina = fcmax(yo), (26) 




FIG. 12: Attraction of the vortex end to the oscillating tip. 
The magnetic tip oscillates with amplitude a = 1.4 along the 
straight line Y — 0.8 parallel to the x axis, starting from 
X = at time t — 0. When the tip approaches the starting 
point from large positive Y, the vortex end shifts to yo ~ 0.11, 
attracted by the tip. At t > the vortex end oscillates along 
x with small, slightly increasing amplitude, moving slowly to 
higher y values. When yo ~ 0.3 is reached the vortex end 
jumps in a few big leaps to its maximum yo ~ 0.73. After 
that it oscillates on a stationary curve. The lower plot shows 
the temporal dependences of xo and yo- All lengths in /im, 
the unit of t is a quarter of the tip period, the parameters are 
as in Fig. 1101 but for simplicity we take f = 1 here. 



aspect ratio is therefore 

max(yo) 

r = 

max(xo) 



<F„ 



F„ 



£>1, 



(27) 



F p the ratio r di- 



and hence max(yo) = F m sin a/k = « / F^ — F^/k. The 



and it is independent of A:. If F„ 
verges, but in this case the vortex displacements are small 
and become less than the vortex radius which is of the 
order of A ab for MFM. For F m and F p of Fig. formula 
(|2T)l yields the aspect ratio r = a/5/3 ~ 1.29, which is 
indeed close to that found in this figure. 

In Fig.rT0]we analyze one more effect that was observed 
experimentally, see Figs. 4c and 4d in Ref. At the ini- 
tial time t — 0, the straight vortex is at x = y — 0. The 
tip oscillates along x with a large amplitude and slowly 
approaches the vortex from large positive y. At a cer- 
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tain time the end of the vortex abruptly jumps to the 
tip and then begins to oscillate with a large amplitude. 
This effect of a sharp onset of the signal is qualitatively 
reproduced by our Fig. [101 A close look to Fig. [TOl shows 
that the large jump of the vortex end is composed of 
several jumps of width increasing nearly exponentially 
in time. These multiple jumps are even better seen in 
Fig. [IT] that shows how the vortex end moves when the 
tip oscillates along y and slowly moves along x start- 
ing far away from the waiting vortex. As compared to 
the corresponding Figs. [5] HI and ITD1 which are described 
by the same parameters and have a vortex-path aspect 
ratio max(j/o)/max(xo)~ 4, in Fig. [TT] the aspect ratio 
max(xo)/max(yo) ~ 1-3 is smaller than even that for the 
isotropic case (« 2.2) since the pinning anisotropy now 
impedes^ the vortex motion in the x direction. 



0.805 




and does not move towards the tip. A more restrictive 
necessary condition for the vortex motion towards the tip 
is that the vortex can oscillate along x. This condition 
yields 



(Y-yo(0)) 2 




m- (28) 
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From numerical calculations we find, see Fig. 1 1 31 that 
there exists a distinct upper threshold for the distance 
(Y — y o (0)) between the vortex end and the tip-oscillation 
line at which the attraction process can occur, and this 
threshold is close to that given by Eq. (|28p If this 
threshold is indeed determined by the pinning forces and 
the dependence of the driving force F on X — xo and 
Y — yo is known, this effect may allow sensitive mea- 
surements of these pinning forces acting on an individual 
vortex in type-II superconductors. 

When the tip oscillates, it generates currents whose ori- 
entation changes in time near the vortex, and the vortex 
motion towards the tip in Fig. [12] as well as the enhanced 
vortex response in the slow scan direction in Figs. [SKI 
are reminiscent of the so-called longitudinal vortex shak- 
ing effect^ In this effect, in essence, a small ac current 
is superimposed perpendicularly to a dc critical current 
that flows in a sample. This leads not only to a periodic 
tilt of vortices but also to their unidirectional drift along 
the direction of the ac current and causes a dc electric 
field along the dc current. In the considered case of the 
oscillating tip the currents flow only near the surface of 
the superconductor, and only the upper, depinned part 
of the vortex "drifts" . 



FIG. 13: Attraction of the vortex end (xq, yo) to the tip oscil- 
lating along x at constant Y as in Fig. [12] but for various dis- 
tances Y — 0.79, . . . 0.85. Plotted is the maximum value j/ m ax 
of yo in each half oscillation vs. time t. For 0.79 < Y < 0.81 
this j/ ma x is slowly increasing and then suddenly jumps to a 
saturation value « 0.73 within about five half oscillations. At 
Y — 0.81 this steep jump occurs only at t = 1500 (after 375 
oscillations). For Y > 0.815, y ma x saturates exponentially in 
f to a small value < 0.114, and thus there is no jump. All 
parameters and units are the same as in Fig. 1121 

In Fig. [12] we reproduce one more experiment described 
in Ref . |l| (in the supplementary material) . In this exper- 
iment the tip oscillates along a straight line at t > and 
does not shift in the perpendicular direction. The vortex 
that waits at some distance from the line of the tip oscil- 
lations at t < 0, at t > begins to move towards the tip. 
Figure [12] shows this attraction process for the isotropic 
case £ = 1 and for the tip-oscillations line parallel to the 
x axis. The initial shift yo(0) of the vortex end along 
y occurs at t < when the tip approaches its starting 
point from large positive Y. This shift occurs if the driv- 
ing force at t = exceeds the appropriate pinning force 
F p = fpX. If the initial distance of the vortex from the 
tip-oscillations line is so large that the driving force is less 
than this pinning force, the vortex end remains pinned 



IV. CONCLUSIONS 

We derive equations that describe the deformation of 
an individual vortex in anisotropic type-II superconduc- 
tors under the influence of the moving tip of a magnetic 
force microscope. These equations take into account the 
driving force generated by the tip, the elastic force caused 
by the vortex deformation, and the pinning force exerted 
by point defects. These equations are valid even at large 
deformations of the vortex, and they properly allow for 
the biaxial anisotropy of the superconductor. From these 
equations, we reproduce the main features of the experi- 
mental data obtained recently^ In particular, we explain 
the enhanced response of the vortex to pulling in the 
slow scan direction as compared to its response in the 
direction of the fast zigzag scan. We demonstrate that 
the vortex position at time t depends on the trajectory 
of the tip at previous times, and it is this property that 
eventually leads to the enhanced vortex response in the 
slow scan direction. We also point out that the enhanced 
anisotropy of pinning in the a-b plane that was observed 
in Ref. [0 is partly caused by the fact that the critical 
force at which the vortex starts to move depends on the 
angle 9 of the vortex tilt and in general does not coincide 
with the pinning force. 
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We note a still unresolved problem. In order to ob- 
tain quantitative agreement of our calculations with the 
experimental data, we have to take a larger value of the 
vortex line tension than the value following from the the- 
oretical estimate. The small line tension ~ e 2 e of a vor- 
tex in an anisotropic bulk superconductor results from 
the almost complete cancellation of the increase of the 
length of a tilted vortex and the decrease of its energy 
per unit length, ei(0) « £o cos 9, with increasing tilt angle 
9Z The existence of the surface at z = and of the tip 
changes the energy ei (9) in the surface layer of depth A 
and, consequently, the almost complete cancellation does 
not occur there. The line tension of a vortex segment 
near the surface may thus be noticeably larger than the 
tension in the bulk. The discrepancy also may be due to 
the too simple expressions for the lateral driving force, 
see Eqs. JT]) and (JU). Since the penetration depth A of the 
driving force should be of the order of A a &, this A is com- 
parable with the experimental values of Z + ho- In this 
situation the correct driving force acting on a curved vor- 
tex at small distances R ~ Z+ho from the tip, is likely to 
be given by formulas more complicated than Eqs. (JTJ) and 
([2]). But Eq. (fTJ) was, in fact, used in the experiment^ for 
the extraction of the lateral driving force, which might 
lead to some overestimation of this force. Thus, a more 
detailed theoretical investigation of the driving force and 
the nonlocal line tension near the surface is needed. 

One more problem that should be studied both the- 
oretically and experimentally is the vortex-motion ran- 
domness that is superimposed on the regular vortex mo- 
tion considered here. This randomness is clearly seen in 
the experimental data of Auslaender et alX It is quite 
possible that apart from point defects and the weak col- 
lective pinning associated with them, in the sample there 
may be strong pinning centers, e.g., the clusters of point 
defects discussed in Ref. [l|, that lead to the observed ran- 
domness. 
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APPENDIX A: FORMULAS FOR THE CRITICAL 
FORCE 

Taking into account formulas ([7]) and ([8]), we obtain 
the following explicit form of Eqs. (TlTj) and (I12p : 



tan?/> 



Wi) 



1/2 



(Al) 
(A2) 



where f p (ipi) = f p (0, tp, ipi) is given by Eq. (JT]), the prime 
means d/dip\, 

N = tan 3 tpi + 1.5 AC tan 2 V>i cos 9 sin 2tp 

+ tan ipi (2r)i cos 2 9 — rj) — 0.5 AC cos 9 sin 2<p, 

D = rji cos 2 6> + 1.5A£tanV'i cos9sin2tp (A3) 

+ tan 2 ?/>i (2t) - rji cos 2 9) - 0.5AC tan 3 ^i cos 9 sin 2tp, 
| £' | AC cos 9 sin 2(p cos 2ip\ + sin 2ipi (r) — cos 2 6r)i ) 

£ A^ cos 9 sin 2ip sin 2ipi+ 2 cos 2 tp\ cos 2 9r]i+ 2 sin 2 tpi rj ' 

and AC = C~ C 1 ? V = J ?( < / 3 ); Vi = vif + tt/2). Equation 
(|A1[) permits one to find the auxiliary angle ipi in terms 
of ip, and then one can calculate f c from Eq. (|A2j) . 
Let us define the parameter 5 by the formula 



25= 



cos 9 



-rji cos ( 



V( 



cos( 



-771 cos 9 



(A4) 

For example, if <p = 0, one has S = cos 2 9/( 2 , while if 
tp = 7r/2, we obtain S = min[cos 2 9( 2 , (cos 2 9( 2 )~ 1 ]. If 
the parameter S is larger than 1/2, there is a one-to-one 
correspondence between ip and At 5 < 1/2 the sit- 
uation changes. 6 - In this case spurious branches of V'i(V') 
appear. The physical branch corresponds to a minimum 
value of f c . For example, in the case 6 < 1/2 one finds 
the following expression for the critical force at tp = 
and ip = or ip = ir/2: 



f c {i> = 0; cp = 0, tt/2) - 2f p0 VS(l-S), 



(A5) 



where f p0 = f£ VC at tp = and f p0 = /p/VC at tp = tt/2. 
On the other hand, at S > 1/2 one has f c (ip = 0;tp = 
0, tt/2) = f p o instead of formula (|A5[) . 
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